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We predict the locations of several multicritical points of the Potts spin glass model on the 
triangular lattice. In particular, continuous multicritical lines, which consist of multicritical 
points, are obtained for two types of two-state Potts (i.e., Ising) spin glasses with two- and 
three-body interactions on the triangular lattice. These results provide us with numerous 
examples to further verify the validity of the conjecture, which has succeeded in deriving 
highly precise locations of multicritical points for several spin glass models. The technique, 
called the direct triangular duality, a variant of the ordinary duality transformation, directly 
relates the triangular lattice with its dual triangular lattice in conjunction with the replica 
method. 

KEYWORDS: spin glass, multicritical point, Potts model, exact solution, triangular lattice, 
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1. Introduction 

Numerous investigations and analyses have been carried out for many years to study 
the properties of finite-dimensional spin glasses. While there have been many studies based 
on mean-field analyses, approximate approaches, and numerical simulations, there are much 
fewer exact or rigorous approaches. 1-3 Among these studies for finite-dimensional spin glasses, 
a reliable approach has been developed to predict the exact location of the multicritical point, 
we call this approach the conjecture. Since its proposal, many predictions for the locations of 
the multicritical points have been given and verified by numerical studies, including the ± J 
Ising model on the square, 4 ' 5 triangular, and hexagonal lattices, 6 ' 7 and various mutually dual 
pairs of hierarchical lattices. 8 

While investigation of the multicritical point is an issue in the phase diagram of spin glass 
models, it has also been pointed out that the location of the multicritical point of the ± J Ising 
model is closely related to an accuracy threshold of reliability of the toric code, one of the 
error-correcting codes used in quantum information. 9 It is therefore an important problem 
to determine the location of the multicritical point, not only for its intrinsic interest, but 
also for its influence in other branches of physics. Most of the predicted results are in good 
agreement with numerical estimations, within their error bars. However, the validity of the 
conjecture is still uncertain, and small deviations have been observed in a few samples. 8 We 
must therefore further verify the validity of the conjecture for many examples and clarify the 
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required conditions for obtaining the exact location of the multicritical point. This is the goal 
of the present work. 

Here, we introduce the Potts spin glass model with two- and three-body interactions. The 
pure, non-random Potts model with two- and three-body interactions satisfies self-duality 
under the direct triangular duality. 10 Because this model has two independent coupling pa- 
rameters, the phase boundary is described as a line in the two-dimensional phase diagram. 
Similarly to this pure Potts model, the Potts spin glass model may have a line of multicritical 
points. Based on these predictions, we investigate numerous examples to check the validity 
of the conjecture. It may also be possible to investigate in detail the required conditions 
for obtaining the exact location of the multicritical point and to find the reason why some 
discrepancies have been observed. 

This paper is organized as follows. We first review the properties of the pure Potts model 
with two- and three-body interactions on the triangular lattice in §2. We introduce, in §3, the 
two types of Potts spin glass models with two- and three-body interactions on the triangular 
lattice. In §4, we discuss the duality for these Potts spin glass models and the possibility of 
predicting their multicritical points. We show that several cases of the Potts spin glass models 
satisfy the known conditions for obtaining the exact location of the multicritical point: the 
2-state Potts (i.e., Ising) spin glass models, and the g-state Potts spin glass models only 
with three-body interactions. We show the results of the presented analyses in §5. The results 
predict multicritical lines, consisting of multicritical points, for the two types of Ising spin glass 
models and multicritical points for the g-state Potts spin glass model only with three-body 
interactions. Finally, we discuss the present work. 

2. Potts model with two- and three-body interactions 

The Potts model with two- and three-body interactions on the triangular lattice is a gener- 
alization of the ordinary Potts model with only two-body interactions on the triangular lattice. 
It was first proposed by Kim and Joseph through the duality and star-triangle transformations 
for the ordinary Potts model on the triangular lattice, which has only two-body interactions 
on each bond. 11 Subsequently, the direct triangular duality was developed by simultaneous use 
of the duality and star-triangle transformations. 10 This direct triangular duality was applied 
to the Potts model with two- and three-body interactions as a general case of the analysis by 
Kim and Joseph. It will be helpful to first review the pure g-state Potts model with two- and 
three-body interactions on the triangular lattice before we consider the random systems. 

We denote the spin variables as fa and the differences as faj = fa — 4>j , which take integer 
values from to q — 1 (mod g). Then, the partition function of the Potts model with two- and 
three-body interactions is written as 




(1) 
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Fig. 1. The Potts model with two- and three-body interactions on a triangular lattice. Each of the 
spin variables interacts with nearest-neighboring sites (two-body interactions). We also consider 
three-body interactions among spins on all up-pointing triangles, shown as shaded. No three-body 
interactions are considered on the down-pointing triangles. 

Here, the product over A runs over up-pointing shaded triangles shown in Fig. 1, and A[- ■ ■ ] 
is the face Boltzmann factor on each up-pointing triangular face: 



where i ^ j runs over the bonds surrounding each up-pointing triangle (in particular, (12), 
(23), and (31) in Fig. 1). if 2 and K% are coupling constants of the two- and three-body 
interactions. Thus, we describe the thermodynamic state as a point on the K2-K3 plane. 
In the usual definition of the partition function, K2 and K% are proportional to the inverse 
temperature (3 = 1/ksT. From this viewpoint, a line is described on the K2-K3 plane by a 
change in temperature. This line is called the thermodynamic line. 12 If the thermodynamic 
line goes across the phase boundary, the Potts model undergoes a phase transition. 

Next, we consider how to determine where the phase transitions occur. The duality trans- 
formation is a useful tool to give the exact location of the transition point. However, an 
ordinary duality, which is carried out by the Fourier transformation, 13 changes the structure 
of the lattice. For example, in our problem, the triangular lattice is transformed into a hexag- 
onal lattice. Therefore, unlike the simple case of the pure Ising model on the square lattice, it 
is necessary to apply another transformation, the star-triangle transformation, to relate the 
hexagonal lattice with another triangular lattice. In this paper, we call the simultaneous use of 
these techniques the direct triangular duality. It is carried out by application of the following 
transformation: 



A« 2 ' K3 [<f> 12 , fe, 0ai] = exp I K 2 Sttij) + K 3 YlSifrj) \ , 



(2) 




{4>ij } 
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x exp | — — (fci0 23 + k 2 (j)31 + h<p 12 ) | , (3) 

where A and A* are the original and dual face Boltzmann factors, respectively. The summation 
with the subscript {<f>ij} runs from {0i2, 023, 03i} = {0,0,0} to {q — l,q — l,q — 1}. The 
quantities ki are spin variables on the dual triangular lattice, and kij represent their differences 
ki — kj. In eq. (3), the exponential terms are produced by the ordinary duality, and Kronecker's 
delta arises from the application of the star-triangle transformation. 

The face Boltzmann factor (2) can be expanded by the identity exp(K6(- ■ ■ )) = l+v5(- • • ), 
where v = ex.p(K) — 1, as follows: 

A K ^ KZ [012, 023, <hi] = 1 + v 2 £ Sfaj) + y 1] Sfaj). (4) 

Here, y = 3v\ + v\ + v 3 (l + v 2 f = exp(K 3 + 3K 2 ) - 3exp(K 2 ) + 2, v 2 = exp(K 2 ) - 1, 
and v 3 = exp(Ks) — 1. By the direct triangular duality (3), the face Boltzmann factor is 
transformed to 

A%' K ° [k 12 , k 23 , k 31 ] = V - I 1 + ^ + - II Kkj) \ ■ (5) 

This equation states that the Potts model with two- and three-body interactions is a self-dual 
model which satisfies 

v 2 = , y = —■ (6 

v y 

The fixed-point conditions are therefore v 2 = v 2 and y = y*. Both of these fixed-point condi- 
tions give a single equation, y = q. This condition y = q gives a line consisting of the fixed 
points. This fixed line is expected to determine a part of the phase boundary of the Potts model 
with two- and three-body interactions. Wu and Zia 12 argued that this fixed line indeed is a 
phase boundary when the couplings satisfy the criteria J 3 + 3J 2 > and J3 + 2 J 2 > 0. These 
criteria mean that coupling constants K 2 and K 3 make effective ferromagnetic interactions; 
i.e., all the spins should be parallel in the ground state. 
When v 2 = 0, the duality relations (6) are reduced to 

«S = -■ (7) 

^3 

In other words, the duality relations (6) map the Potts model with two- and three-body 
interactions (v 2 / 0, v 3 ) into another Potts model, that is, the dual model, with two- and three- 
body interactions (v 2 / O,^). If we apply the direct-triangular duality to the Potts model 
with only three-body interactions (v 2 = 0,^3), it is transformed into the dual Potts model 
with only three-body interactions (v 2 = 0,w|). It is obvious that the fixed-point condition for 
the Potts model only with three-body interactions is v 3 = q. 
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3. Potts spin glass model 

In the previous section, we introduced the Potts model with two- and three-body inter- 
actions. Here, we add quenched randomness to the two-body and three-body interactions in 
this Potts model. 

The Hamiltonian of the Potts spin glass model on the triangular lattice is defined as 



i+3 i¥=3 



m 



i j j 



and its face Boltzmann factor is 



A* 2,Ki [{4>ij + kj}; {4>ij + rriij}] = exp ! K 2 ^ Sfaj + kj) + K 3 <5(<^ + m 



(8) 



(9) 



i¥=3 



&3 



where faj are the same symbols as in the previous section. The quantities lij are the random 
variables for the two-body interactions, and rriij are those for the three-body interactions. 
These random variables follow the distribution functions 



l-(g-l)P2 I 







P{{m i3 }) = { 



p 2 otherwise 
= p 2 ex.p{K p2 5(li j )} 

l-(q 3 -l)P3 {my} = {0,0,0} 
j>3 otherwise 



(10) 



(11) 



= P3exp{if p3 JJ^m^)}, 

where exp(K p2 ) = {I — (q — I)p2}/P2 and exp(K p3 ) = {1 - (q 3 - l)j> 3 }/;P3- The brackets {•} 
express a bond set, like (12), (23), and (31) in Fig. 1, on the elementary up-pointing triangles. 

The Hamiltonian of our Potts spin glass model has the property of gauge invariance, 
similarly to the ± J Ising model. The gauge transformation is defined as follows: 



i'.i 



<f>ij ~ (Si ~ Sj) 
kj + ( s i — Sj) 



mi 



> rriij ~i~ (^i Sj)i 

where Sj takes an integer value from to q— 1. The Hamiltonian is invariant under this gauge 
transformation. Considering the form of the distribution functions and the property of this 
gauge invariance, we can calculate the exact internal energy by the well-known procedure 14 ' 15 
in a special subspace, defined by K 2 = K p2 andK 3 = K p3 . This subspace is usually called the 
Nishimori line. Because two independent variables exist in our model, the special subspace 
should be called the Nishimori "surface" in a four-dimensional space (K 2 , K p2 , K3, K p s). It is 
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expected that there are intersections between the Nishimori surface and the phase boundary, 
that is, a multicritical "line" , which consists of the multicritical points. 

It is also possible to define another type of Potts spin glass with two- and three-body 
interactions by using the same random variables for the three-body interactions as those 
for the two-body interactions {kj} instead of the random variables {rriij} in eq. (8). The 
Hamiltonian is then written as 



H = ^ H A [(f) 12 , 023, <fol] = ^ < 



- j 2 $(<pij + kj) - h n si&j + Uj) > • (12) 

In this case, the following distribution function is required if we construct the Nishimori 
surface: 



P({Uj}) = Cexp< 



where C is the normalization factor. If we define the gauge transformation as <pij — > <pij — (sj — 
Sj) and lij — > + (si — Sj), the Hamiltonian is invariant under this gauge transformation. The 
Nishimori surface K 2 = K p2 and K% = K p s is also obtained. We call the former Potts spin 
glass model type I, and the latter one type II. The two slightly different spin glass models are 
equal to each other for the cases only with three-body interactions and only with two-body 
interactions. The following analysis applies to both the type I and type II models, but only 
the type I case is shown. 

4. Duality and conjecture 

As used in the ± J Ising model on the triangular lattice, 7 it is straightforward to develop the 
direct triangular duality for the replicated system. It is possible to apply the direct triangular 
duality to the replicated system as the following replicated duality relation: 



,a=l 



• 27r V^/, 

x exp ^ i — y [k 



a i a 



1 V23 + *2 V31 



Q 



+ k 



3*12) 



(14) 



a=l 



where a is the replica index running from 1 to n, and A n [{</>i 2 }, {023 }> {^31}] i s the 
configurational- averaged and replicated face Boltzmann factor of the type-I Potts spin glass, 
which is given by 

^[{^12}, 

n 

{kj,mij} a=l 

Following the existing successful cases, 5-7 we need two conditions to derive the multicritical 
points: 
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(i) On the Nishimori line or, in our model, surface, the n = 1 and 2 replicated systems are 
self-dual. 

We concentrate on the system on the Nishimori line, where it is believed that the multi- 
critical point lies. 14 ' 15 

(ii) Transition points for replicated systems are determined only by a single equation: 

1*[{0},{0},{0}] =I n [{0},{0},{0}]. (16) 

This kind of equation is satisfied in pure systems, such as the Ising model and the (/-state 
Potts model on the square, triangular, and hexagonal lattices. Indeed it is possible to 
obtain the fixed-point condition y = q for the case of the pure Potts model with two- and 
three-body interactions by substituting eqs. (4) and (5) into eq. (16). It is expected that 
eq. (16) represents the location of the singularity of the replicated system. The validity 
of this hypothesis is still unproven, however. 

We now examine the self-duality for the n = 1 and 2 replicated systems of the type-I 
and type-II Potts spin glass models, as discussed below. We find that the (/-state Potts spin 
glass model with only three-body interactions (v2 = 0, p2 = 1/q) satisfies the self-duality for 
general values of q when the replica numbers are n = 1 and 2. On the other hand, we find that 
only the two-state Potts spin glass models with two- and three-body interactions, for both of 
the type-I and type-II cases, satisfy the self-duality when n = 1 and 2. Neither the type-I nor 
type-II Potts spin glass model with q > 3 states satisfies the self-duality when n = 2. 

Therefore, using the conjecture, it would be possible to obtain the multicritical points 
of the g-state Potts spin glass model only with three-body interactions and both the type-I 
and type-II two-state Potts spin glass models with two- and three-body interactions. On the 
other hand, the conjecture may give incorrect predictions for the Potts spin glass model with 
two- and three-body interactions with q > 3 states. The duality is not always a good tool, 
even for the pure model, to predict the transition point for most systems with many states, 
such as the (/-state clock model (q > 5), which is not self-dual; on the other hand, the (/-state 
clock model is self-dual when q < 4. This property of the duality may be the reason why the 
above unfavorable situation for the replicated models appears when the number of states and 
replicas is higher than a threshold. 

The rest of this section is devoted to a detailed explanation of the above-mentioned facts 
about the self-duality for the replicated models. We consider the self-duality of the type-I 
Potts spin glass model with two- and three-body interactions. 

In the case of the n = 1 replicated Potts spin glass model, the face Boltzmann factor is 
rewritten, from eq. (15), as 

M [012,023,03l] 
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= P 3 2P3(Q 3 + 2v 3 )(q + 2v 2 f 



1 + ^i>«' 




(17) 

If we look at eq. (4) , we find that this is equivalent to the face Boltzmann factor of the Potts 
model with two- and three-body interactions with the effective couplings v~ 2 = v|/(<7 + 2v 2 ), 
and v~ 3 = f|/(<7 3 + 2v 3 ). It is therefore easy to obtain the following fixed-point condition, 
similarly to the case when we obtained the fixed-point condition (6): 

3vj + 4 + v 3 (l+v 2 ) 3 = q. (18) 

In the case of the n = 2 replicated Potts spin glass model, the face Boltzmann factor is 
rewritten, from eq. (15), as: 

= PsP 3 2 WU + ^ + vlj^ 8{4%) + v\b^\ - 4) + v\ f[ 5(^) | 

ijtj I Q=l a=l ) 

g 3 + 3V3 +viY,u m 3 ) + 4 n 6 ^h - 4) + 4 n n w \ ■ ^ 

This face Boltzmann factor does not generally show the self-dual property. The reason is that, 
in eq. (19), the first product over i / j produces the products of Kronecker's delta symbols, 
such as 5{4>\ 2 )&{4% 3 )&{ ( t ) 3 i — 4> 31 ), and its cyclic permutations of the indices. The dual terms of 
these products are given by the application of eq. (14). For example, 5(4 > \ 2 )3(4 ) ' 23 )<5(4 > 3 i ~~ ^ii) 
is transformed as follows: 

^12)^23)^31 " <f>li) - - q 5(k\ 2 - (20) 

However, this type of interaction does not exist in eq. (19). Therefore, the n = 2 replicated 
Potts spin glass model with two- and three-body interactions is not generally self-dual. We thus 
conclude that the condition (i) is not always satisfied for the Potts spin glass model with two- 
and three-body interactions. We find also that the (/-state Potts spin glass model only with two- 
body interactions does not satisfy the condition (i) either, because this problem is produced 
by the first product over i ^ j in eq. (19), which is obtained from terms representing two-body 
interactions in eq. (15). On the other hand, when we consider the case v 2 = (p 2 = 1/q), 
only with three-body interactions, such troublesome terms as 5{4>l 2 )8(4 > 23 )d( ( t ) 3 i — ^31) vanish. 
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More precisely, the face Boltzmann factor (19) is reduced to 

M [{0M,{</>23},{«] 



' , 2 



P3(g 3 + 3V 3 ) < 



' ' <l'> + 3 3« 3 ^n^v ;j, 

= 1 ly^J 



Though this face Boltzmann factor is still a little complicated, the terms involving Kronecker's 
delta symbols, such as Yli^j S(4>fj), and the constant 1 are changed in the following manner 
by the application of the direct triangular duality (14): 

2 

1 - q 2 H5(k«) 8(4%) - «(*g) (13 + a) 

a=l 

2 



Q=l ^ 

From these transformations, the explicit form of the dual face Boltzmann factor is 
^2 [{^12J"> {^23}) {^31}] 

-¥(i4£n«*§>+ s2 ^ 2i n «* - 4) + £ n n 1 ■ 



(22) 



Therefore, the n = 2-replicated g-state Potts spin glass only with three-body interactions is 
self-dual. Thus, it is confirmed that the condition (i) is satisfied for the Potts spin glass model 
only with three-body interactions. The fixed-point condition is expressed as a single equation: 



= q l - (23) 

If we restrict ourselves to the case q = 2, which is equivalent to the Ising system, we find 
a possibility that the system is self-dual. Such terms as S(k\ 2 — k 23 ) are rewritten as follows, 
only in the q = 2 case: 

S(kl 2 ~ k 2 23 ) = 2S(kl 2 )5(k 2 23 ) - 8(k\ 2 ) - 5(k 2 23 ) + 1. (24) 

The terms on the right-hand side can be obtained from eq. (19). Terms consisting of the 
products of the Kronecker's delta symbols not appearing in eq. (19) are obtained by the 
direct triangular duality (14). The n = 2-replicated 2-state Potts spin glass is therefore self- 
dual. In Appendix A, this self-duality is shown more explicitly than here, and the fixed-point 
condition for the two-state Potts spin glass model is derived. The condition (i) is also satisfied 
by the two-state Potts spin glass model. 
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We conclude that it is possible to predict two series of the locations of the multicritical 
points: the multicritical points of the general (/-state Potts spin glass model only with three- 
body interactions and the multicritical lines of both the type-I and type-II two-state Potts 
spin glass models with two- and three-body interactions. 

5. Multicritical points and lines 

Our next tasks are to confirm that the condition (ii) is satisfied and to predict the multi- 
critical points and lines. It is possible to write the following explicit form of eq. (16) for the 
g-state Potts spin glass model, as shown in Appendix B: 

P!P3 {q-l + (l + v 2 ) n+1 f {q 3 - 1 + (1 + « 3 ) n+1 } 3 



P 2 P3 



q\q - lfF n+1 + q 2 {q - l)G n+1 + £ (a^™ +1 + b t G™ +l ) \ . (25) 



j=0 



where 



F = ±-{(q + v 2 f -vl) Fi = F + v 3 {l+v 2 y 

G = ^{{q + v 2 f + (q-l)vl} Gi = G + v 3 (l+v 2 )\ 
and Oj and bi are the coefficients of Fi and Gi, which are evaluated as follows: 

a = (q - l) 3 + (q - l)(q - 2) b = (q-l)(q-2) 

ai = 3(g-l)(g-2) h = 3(q - 1) 

a 2 = 3(g-l) b 2 = 

a 3 = 63 = 1. 

Detailed calculations are shown in Appendix B. Substituting n = 1 and n = 2 with q = 2 into 
eq. (25), we obtain the same equations as the fixed-point conditions for the n = 1- and n = 2- 
replicated Potts spin glasses with two states. When v 2 = (p 2 = l/q), eq. (25) also gives the 
fixed-point conditions (18) and (23). It is confirmed that the condition (ii) is also satisfied 
for the (/-state Potts spin glass model only with three-body interactions and the two-state 
Potts spin glass model. In accordance with the replica method, we consider the n — > limit 
of eq. (25), which is expected to determine the singularity of the Potts spin glass model. We 
obtain the following equation, which gives predictions of the multicritical points for various 
combinations p 2 and p 3 , by taking the leading term of n in eq. (25): 

3 

\q - l) 2 F \nF + q 2 (q- l)G In G + (a^ In F, t + bid In Gi) 



■ 1 </ — 1 ; r 111 f r (/ I, (/ — 1 ,ni 111 r 

-(1 + v 3 )(q + v 2 f ln(l + v 3 ) - 3{q + v 2 f(l + v 2 ){q 3 + v 3 ) ln(l + v 2 ) 



(q 3 + v 3 )(q + v 2 flnq. (26) 
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Fig. 2. The two solid lines connecting (p 2 ,p 3 ) = (1/2,0.0396336) and {p 2 ,p 3 ) = (0.164194, 1/2 3 ) are 
the multicritical lines. The two points T and D are the multicritical points of the Potts spin glass 
model only with three-body interactions (K 2 = 0, p 2 = l/q) and only with two-body interactions 
(K 3 = 0, p 3 = 1/g 3 ), respectively. 

The solutions to this equation will describe a multicritical line on the {p 2 ,p 3 ) or (K p2 ,K p3 ) 
plane. The above result is for the type-I Potts spin glass. We also derive the following equation 
of the multicritical line for the type-II Potts spin glass: 

(q - l)FlnF + G 3 lnG 3 - 3u 3 (l + v 2 fln(l + v 2 ) - 3{q + v 2 ) 2 (l + v 2 ) ln(l + v 2 ) 

-(1 + v 3 )(l + v 2 f ln(l + v 3 ) = {(q - 1)F + G 3 } \nq. (27) 

When we consider the case of v 2 = (p 2 = eq. (26) gives the following equation, which 
determines the location of the multicritical points of the g-state Potts spin glass model only 
with the three-body interactions: 

{1 - q 2 (q - l)p 3 } In {1 - q\q - l)p 3 } - {l - (q 3 - l)p 3 } In {l - (q 3 - l)p 3 ) 

-(q 3 -l)p 3 ln P3 = {l-2q\q-l)p 3 }lnq. (28) 

The solutions are, for example, p 3c = 0.0396336 (q = 2), p 3c = 0.0139431 (q = 3), p 3c = 
0.00637851 {q = 4), and p 3c = 0.00342033 {q = 5). 

Both lines described by eqs. (26) and (27) are shown in Fig. 2. This is a reasonable result 
because both lines smoothly connect the multicritical point of the 2-state Potts spin glass 
model only with three-body interactions p 3c = 0.03969336 with that of the 2-state Potts 
spin glass model only with two-body interactions p 2c = 0.164194, which is equivalent to the 
Ising model on the triangular lattice. 7 Here, p 2 expresses the density of the antiferromagnetic 
interactions. 
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6. Discussion 

We introduced two Potts spin glass models, type-I and type-II, on the triangular lattice. 
Similarly to the ± J Ising model, it is possible to find a special subspace of the phase diagram 
where we can calculate the exact internal energy, known as the Nishimori line. The special 
subspace should be a surface, because the Potts spin glass models have two coupling param- 
eters. It is expected that, considering the properties of the Nishimori line, the multicritical 
points lie on a surface. It is also expected that the multicritical points describe a line, be- 
cause the multicritical points are located at the intersections between the phase boundary and 
this surface. Applying the direct triangular duality in conjunction with the replica method 
to the two types of Potts spin glass models, we argued for the self-duality of the replicated 
systems of these models. Taking account of the existing successful examples of the conjecture 
for many spin glass models, we expect that the condition to obtain the exact location of the 
multicritical points is the self-duality of the replicated models with n = 1 and n = 2 repli- 
cas. We showed the self-duality in only the cases of the g-state Potts spin glass model with 
only three-body interactions and two types of 2-state Potts spin glass models with two- and 
three-body interactions. Thus, our results possibly give the exact location of the multicritical 
points of the general Potts spin glass model with only three-body interactions with q states 
and the multicritical lines of two slightly different Potts spin glass models for the 2-state case. 
However, the hypothesis that a single equation will determine the singularity of the spin glass 
models is not completely verified. Therefore, we have to compare our results with the results 
obtained by other techniques. 

On the other hand, neither of the n = 2 replicated Potts spin glass models with two- and 
three-body interactions with q > 3 states satisfy the self-duality, which may result in incorrect 
predictions for this case. However, it is not certain which conditions are truly necessary. It 
would also be worth investigating the location of the multicritical lines of the Potts spin glass 
models with q > 3 states to verify the necessity of the self-duality for the n = 2-replicated 
Potts spin glass model. 

For the above reasons, it is expected that our multicritical lines with 2 states will be 
helpful to further investigate the validity of the conjecture. Moreover, an investigation of the 
phase transition of the Potts spin glass model with q > 3 states by other approaches should 
elucidate the precise conditions necessary to obtain the multicritical points. 

7. Acknowledgements 

The author would like to thank Prof. H. Nishimori of the Tokyo Institute of Technology for 
useful discussions and suggestions on the manuscript. This work was partially supported by 
the 21st Century COE Program 'Nanometer-Scale Quantum Physics' at the Tokyo Institute of 
Technology, and by Grant-in- Aid for Scientific Research on the Priority Area "Deepening and 
Expansion of Statistical Mechanical Informatics" from the Ministry of Education, Culture, 



12/18 



J. Phys. Soc. Jpn. 

Sports, Science, and Technology of Japan and by CREST, JST. 



Full Paper 



(A-l) 



Appendix A: Duality of the n = 2-replicated Potts spin glass model with 2 states 

It is a little complicated to write the expansion of eq. (19). Therefore, we directly calculate 
the face Boltzmann factors considering symmetry in eq. (19). Because of the constraint (j>f 2 + 
023 + 4>si = w ith mod 2, the following four types of face Boltzmann factors from eq. (19) 
characterize the replicated Potts spin glass model: 

A) = Ao[0, 0,0; 0,0,0] 

= P&3 (q + 3u 2 + 3vj + v%) 3 (q 3 + 3v 3 + 3vj + vf) 

i 2 [0, 0,0; 0,1,1] = A 2 [0,0,0; 1,0,1] = A 2 [0, 0, 0; 1, 1, 0] 

i 2 [0, 1, 1; 0, 0, 0] = A 2 [l, 0, 1; 0, 0, 0] = A 2 [l, 1, 0; 0, 0, 0] 

p\pz {q + 3v 2 + Zvl + vfj (q + 3v 2 + vl) 2 (q 3 + 3v 3 + vfj 

i 2 [0, 1,1; 0,1,1] =i 2 [l, 0,1; 1,0,1] = A 2 [l, 1,0; 1,1,0] 

plP3 (q + 3f 2 + 3v 2 2 + vf) (q + 3v 2 + vf) 2 (q 3 + 3v 3 ) 

i 2 [0, 1,1; 1,0,1] =i 2 [0, 1,1; 1,1,0] =i 2 [l, 0,1; 0,1,1] = A 2 [l, 0,1; 1,1,0] 
i 2 [l, 1,0; 0,1,1] =i 2 [l, 1,0; 1,0,1] 



A = 



A = 



As = 



(A-2) 



(A-3) 



= Pip 3 (g + 3f 2 + ^) 3 (g 3 + 3t;3). 



(A-4) 



Using these face Boltzmann factors, it is possible to write its dual face Boltzmann factor using 
eq. (14), as follows: 

A 2 [{kf 2 }, {^23}, {^31}] 
2 

A + Aj^J^exp (iTTfcg) 



a=1 



(A-5) 



+A2 ^ exp {in^ + 4)} + A 3 exp {i7r(fcj fc + k 2 ki )} 

ijtj ijk 

This face Boltzmann factor has the same symmetry as the original one AK012/) {^3}; {^31}] 
for example, A* 2 [0, 0, 0; 0, 1, 1] = A* 2 [0, 0, 0; 1, 0, 1] = A?[0, 0, 0; 1, 1, 0] = A*[0, 1, 1; 0, 0, 0] = 
A;[l, 0, 1; 0, 0, 0] = A 2 [l, 1, 0; 0, 0, 0]. Therefore, the following duality relations are obtained: 

An 



A*, 



- (An + 6 A + 3A2 + 6 A) 

- (An + 2A - A2 - 2A) 



^2 = 4 (A - 2A + 3A - 2A 3 ) 



A = 



(An - 2A X - 1A 2 + 2 A3) • 



(A-6) 
(A-7) 
(A-8) 
(A-9) 
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Setting u>i = A^/Aq, the above equations are rewritten as 

* 1 



Ui + - 



1 + 6cji + 3^2 + 6c<J3 



(A-10) 



It is therefore possible to obtain the following fixed condition, where = uj* for any i: 

1 + 6wi + 3^2 + 6o; 3 = 4. (A-ll) 

This condition can again be derived by A = Aq, that is, eq. (16). It is also straightforward 
to show that the general cases of q states are not self-dual, similarly to the above discussion, 
or by direct manipulation of eq. (14). 

Appendix B: Evaluation of eq.(25) 

We will show the evaluation of eq. (25) in this section. The quantity of eq. (25) is obtained 
by substituting {(f>fj} = {0} into eq. (15) and by substituting {kfj} = {0} into eq. (14) under 
the Nishimori surface conditions K<i = K p 2, and K% = K p %. 

First, the left-hand side of eq. (25) is obtained from eq. (15), as follows: 

i n [{0},{0},{0}] 



pIpsJ2 E a a 3 [{%>; Wl >< II • * 3 [{U;{"H>>] 

{hj}{mij} a=1 



: pIps e e n exp {( n+i ) K2<5 (^)} exp< 

{hj} {mij} i+3 



(n + l)K 3 l[5(m i3 ) 

i+3 



^ j 1 + (e (n + l)K 2 



.{hj} 



J2 \ l + (e^ K3 -l)l[S( 



m. 



: p\pz {? - 1 + (1 + v 2 ) n+1 Y {q 3 - 1 + (1 + v 3 ) n+1 } . 



(B-l) 



Next, we derive the expression on the right-hand side of eq. (25). We substitute {k?j} = 
{0}, that is {kf } = {kf}, into eq. (14) under the conditions of the Nishimori surface. Then, 
the exponential term in eq. (14) vanishes because of the constraints by Kronecker's delta 
5((f)i2 + 023 + ^31 ) f° r eacn replica. It is therefore necessary to evaluate only the following 
quantity: 

i*[{0},{0},{0}] 



P2P3 



J2 E ^a'* 3 [{U;{^}] 



{kj}{mij} 



n 

a=l 



E ^2 + 023 + <f>si)^' K3 im + kjY, m + mij}] 

,{4>ij} 



(B-2) 
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n 



a=l 



E ^ + $3 + ^l)< 2 ' K3 m + kjh {<P% + my}] 



n 



E ^ + 023 + <t%l) Ili 1 + ^<K<fe + ^)H1 + «3 II WiJ + ^ 



= F + C({l ij };{m ij }), 
where the two quantities F and C({hj} : {my}) are defined as follows: 
F = - {{q + v 2 f - vl) 



(B-3) 



(B-4) 



C({kj}; {rriij}) = v 3 5(mi 2 + m 23 + m 31 ) + v 2 S(lij - rriij)} + v|<5(Zi 2 + l 23 + foi)- 

(B-5) 

Substituting eq. (B-3) into eq. (B-2) yields 
i*[{0},{0},{0}] 

= pIpsv 3 E II i 1 + V* 6 ^)} {F + C({lijh {0})}" 

{hj } i+3 

+PIP3 E E II + ^fe)} ^ F + C ^ W)F • (B'6) 

We divide both sets {hj} and {my} into two subsets P and P to sum over {hj} and {my}. 
Here, P is a subset which satisfies the constraints X\ 2 + x 23 + x 3 \ 7^ 0(mod q), X{j G P. The 
reason why we consider such division of sets {hj} and {my} is that the quantity C({Zy}; {my}) 
depends on the values of l\ 2 + l 23 + /31 and mi2 + m 23 + m3i. Dividing sets {hj} and {my} 
and expanding the terms with power n yields 



^[{0},{0},{0}] 



P 2 P3V3 



E E 

{Jy}gP r=0 



n 
r 



F"-^nexp{(r + l)^ 2 <5(Zy)} 



+ E E 



n 
r 



(F + « 2 3 )^ r vl [] exp {(r + l)K 2 5(h,)} 

i¥=3 



P 2 P3 



E E F n Y[e W {K 2 5(h,)} 

{hj}eP {mij}eP i+j 

+ E E n ex p{^<5(%)} 

{kj}£P {mij}&P i+j 
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+ £ £ £ 

n 

+ £ £ £ 

{kj}eP {m i:j }eP r=0 



n 
r 



n 
r 



F n ~ r vl Yl exp {K 2 6(Uj) + rK 2 S{Uj ~ mj)} 

(F + vl) n ~ r v r 3 Yl exp {Krftkj ) + r J PC 2 <5(Z ij - m^)} 



(B-7) 



It is possible to evaluate all terms except for the last two terms in this equation as follows: 



P 2 P3 



v 3 + V2Y + biG?(l + v 2 y} + q\q - lfF n+l + q\q - l)G 



n+1 



(B-8) 



(B-9) 



q" 
where 

a = (q - l) 3 + (q - l)(q - 2) 6 = (g - 1)(? - 2) 

ai = 3(g-l)(g-2) h = 3(q - 1) 

a 2 = 3(g-l) 6 2 = 

a 3 = 63 = 1. 

Here, we used the following convenient identities: 
^[]exp{i?%)} = (q + v R f 

{hj} i¥=j 

HexpimUj)} = ^-{(q + v R f + (q-l)v R }=G 

n^pw^)} = (i-^){(9+^) 3 -4}=(9-i)^- (B-i2) 

On the other hand, the last two terms of eq. (B-7) are evaluated as follows, considering the 
relationship of the two subsets P and P: 



(B-10) 
(B-ll) 



P2P3 



E E E 

.{lij}eP{mij}ePr=0 



n 
r 



+ E E E 

{k 3 }eP{ mi] }epr=o 



n 
r 



p\P2, 



E E E 



n 
r 



(F + t>l) n " r ^ J] exp {K 2 <5(^) + rK 2 8{Uj - m^)} 
j F"- r «5 J] ex P {^2*(iij) + rK 2 bil id - rmj)} 



+ £ £ £ 

{hj}&P {n Hj }&P r=0 



n 
r 



(B-13) 



16/18 



J. Phys. Soc. Jpn. Full Paper 

We use another expression for limitation of the summation over {hj} € P and {m^} £ P: 

Yl = Y + hs + hi) 

{kj}&P {hj} 

9-1 



- y~\ yz ex p 1^—^12 + + ki) \ ■ 



Eq. (B-13) is then evaluated as follows 

q— 1 n 



i EE 1 p m - r «5{«(9 + «a + t4 r) )«(7,) + ^} { 



r,=0 r=0 
(J— 1 9—1 n 



4eee " {(f+^-f-}- 



r;=0r m =0 r=0 



|g 2 <5(r;)5(r m ) + qv 2 5(Ti) + qv^ 5{ji + r m )w 2 + 4 r) } 



(B-14) 



where = exp^i^) — 1- It is straightforward to expand this equation and obtain the 
following expression: 



P2P3 



3 3 
J=0 i=0 



(B-15) 



Adding eq. (B-8) to eq. (B-15), we obtain the following result for A* n [{0}, {0}, {0}], which is 
the right-hand side of eq. (25): 



i;[{o},{o>,{o}] 



q\q - \fF n+1 + q\q - l)G n+1 + £ a^ n+1 + £ b t G^ +1 



i=0 



i=0 



(B-16) 
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